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ABSTRACT 

Monitoring  surface  currents  by  coastal  high-frequency  radars  (HFRs)  is  a  cost-effective  observational  tech¬ 
nique  with  good  prospects  for  further  development.  An  important  issue  in  improving  the  efficiency  of  HFR 
systems  is  the  optimization  of  radar  positions  on  the  coastline.  Besides  being  constrained  by  environmental  and 
logistic  factors,  such  optimization  has  to  account  for  prior  knowledge  of  local  circulation  and  the  target  quantities 
(such  as  transports  through  certain  key  sections)  with  respect  to  which  the  radar  positions  are  to  be  optimized. 

In  the  proposed  methodology,  prior  information  of  the  regional  circulation  is  specified  by  the  solution  of  the 
4D  variational  assimilation  problem,  where  the  available  climatological  data  in  the  Bering  Strait  (BS)  region 
are  synthesized  with  dynamical  constraints  of  a  numerical  model.  The  optimal  HFR  placement  problem  is 
solved  by  maximizing  the  reduction  of  a  posteriori  error  in  the  mass,  heat,  and  salt  (MHS)  transports  through 
the  target  sections  in  the  region.  It  is  shown  that  the  MHS  transports  into  the  Arctic  and  their  redistribution 
within  the  Chukchi  Sea  are  best  monitored  by  placing  HFRs  at  Cape  Prince  of  Wales  and  on  Little  Diomede 
Island.  Another  equally  efficient  configuration  involves  placement  of  the  second  radar  at  Sinuk  (western 
Alaska)  in  place  of  Diomede.  Computations  show  that  1)  optimization  of  the  HFR  deployment  yields  a  sig¬ 
nificant  (1.3-3  times)  reduction  of  the  transport  errors  compared  to  nonop timal  positioning  of  the  radars  and  2) 
error  reduction  provided  by  two  HFRs  is  an  order  of  magnitude  better  than  the  one  obtained  from  three 
moorings  permanently  maintained  in  the  region  for  the  last  5  yr.  This  result  shows  a  significant  advantage  of  BS 
monitoring  by  HFRs  compared  to  the  more  traditional  technique  of  in  situ  moored  observations.  The  obtained 
results  are  validated  by  an  extensive  set  of  observing  system  simulation  experiments. 


1.  Introduction 

In  recent  years,  the  technology  of  monitoring  coastal  cir¬ 
culations  by  HF  radars  has  improved  considerably  (Kohut 
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and  Glenn  2003;  Kohut  et  al.  2012)  and  now  routinely 
provides  real-time  maps  of  surface  currents  up  to  250  km 
offshore  (Harlan  et  al.  2010;  Barrick  et  al.  2012;  Abascal 
et  al.  2012).  Because  of  limitations  in  the  number  of  con¬ 
tinuously  maintained  high-frequency  radars  (HFRs),  it  is 
desirable  to  optimize  their  locations  in  such  a  way  that  HFR 
readings  could  be  easily  converted  into  robust  estimates  of 
the  key  circulation  parameters  [target  quantities  (TQs)]. 


©  2015  American  Meteorological  Society 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

FEB  2015 


2.  REPORT  TYPE 


3.  DATES  COVERED 


00-00-2015  to  00-00-2015 


4.  TITLE  AND  SUBTITLE 


Optimization  of  the  High-Frequency  Radar  Sites  in  the  Bering  Strait 
Region 


6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory, Stennis  Space  Center, MS, 39556 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Monitoring  surface  currents  by  coastal  high-frequency  radars  (HFRs)  is  a  cost-effective  observational 
technique  with  good  prospects  for  further  development.  An  important  issue  in  improving  the  efficiency  of 
HFR  systems  is  the  optimization  of  radar  positions  on  the  coastline.  Besides  being  constrained  by 
environmental  and  logistic  factors,  such  optimization  has  to  account  for  prior  knowledge  of  local 
circulation  and  the  target  quantities  (such  as  transports  through  certain  key  sections)  with  respect  to  which 
the  radar  positions  are  to  be  optimized.  In  the  proposed  methodology,  prior  information  of  the  regional 
circulation  is  specified  by  the  solution  of  the  4D  variational  assimilation  problem,  where  the  available 
climatological  data  in  the  Bering  Strait  (BS)  region  are  synthesized  with  dynamical  constraints  of  a 
numerical  model.  The  optimal  HFR  placement  problem  is  solved  by  maximizing  the  reduction  of  a 
posteriori  error  in  the  mass,  heat,  and  salt  (MHS)  transports  through  the  target  sections  in  the  region.  It  is 
shown  that  the  MHS  transports  into  the  Arctic  and  their  redistribution  within  the  Chukchi  Sea  are  best 
monitored  by  placing  HFRs  at  Cape  Prince  of  Wales  and  on  Little  Diomede  Island.  Another  equally 
efficient  configuration  involves  placement  of  the  second  radar  at  Sinuk  (western  Alaska)  in  place  of 
Diomede.  Computations  show  that  1)  optimization  of  the  HFR  deployment  yields  a  significant  (1.3???3 
times)  reduction  of  the  transport  errors  compared  to  nonoptimal  positioning  of  the  radars  and  2)  error 
reduction  provided  by  two  HFRs  is  an  order  of  magnitude  better  than  the  one  obtained  from  three 
moorings  permanently  maintained  in  the  region  for  the  last  5  yr.  This  result  shows  a  significant  advantage 
of  BS  monitoring  by  HFRs  compared  to  the  more  traditional  technique  of  in  situ  moored  observations.  The 
obtained  results  are  validated  by  an  extensive  set  of  observing  system  simulation  experiments. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

13 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


298 


JOURNAL  OF  ATMOSPHERIC  AND  OCEANIC  TECHNOLOGY 


Volume  32 


The  theory  of  optimal  observations  has  been  consid¬ 
ered  by  many  authors  (e.g.,  Baker  and  Daley  2000;  Kohl 
and  Stammer  2004)  and  was  used  extensively  in  dynam¬ 
ical  sensitivity  studies  (e.g.,  Marotzke  et  al.  1999;  Losch 
and  Heimbach  2007;  Veneziani  et  al.  2009;  Zhang  et  al. 
2011)  and  optimization  of  in  situ  observations  (Panteleev 
et  al.  2009).  The  key  idea  is  to  explore  the  impact  that 
small  variations  in  the  likelihood  (cost)  function  attrib¬ 
uted  to  new  observation  sites  have  on  the  a  posteriori 
error  variances  of  the  optimized  values  of  the  TQs  (e.g., 
velocity  values,  in  particular  locations  or  the  mass,  heat, 
and  salt  (MHS)  transports  through  particular  sections). 
Observation  sites  providing  the  largest  impact  on  the  TQs 
are  assumed  to  be  optimal  in  the  sense  that  they  provide 
the  most  information  on  the  variability  of  TQs,  and 
therefore  most  effectively  reduce  the  a  posteriori  error 
variance  after  the  corresponding  data  have  been  assimi¬ 
lated  into  a  numerical  model. 

In  the  present  study,  we  explore  application  of  the  four¬ 
dimensional  variational  data  assimilation  (4DVar)  sensi¬ 
tivity  technique  to  the  problem  of  the  optimal  placement 
of  HFRs.  Performance  of  the  method  is  assessed  in  the 
Bering  Strait  (BS)  region,  which  is  important  from  both 
academic  and  practical  viewpoints.  First,  the  anticipated 
growth  of  pan-Arctic  transportation  requires  accurate 
assessment  and  forecasts  of  local  circulation.  From  the 
academic  point  of  view,  the  BS  and  the  Chukchi  Sea  (CS) 
comprise  the  “hot  spot”  controlling  the  exchange  of  ice 
and  seawater  between  the  Pacific  and  Arctic  Oceans.  The 
impact  of  BS  fluxes  on  the  thermohaline  structure  of  the 
Arctic  Ocean  (AO)  depends  on  the  redistribution  of 
the  Pacific  heat  and  freshwater  anomalies  within  the  CS 
basin  (Fig.  1).  Both  observational  (Woodgate  et  al.  2005; 
Weingartner  et  al.  2005;  Luchin  and  Panteleev  2014) 
and  modeling  (Winsor  and  Chapman  2004;  Spall  2007; 
Watanabe  and  Hasumi  2009)  studies  suggest  that  Pacific 
waters  enter  the  open  AO  through  two  pathways  (Fig.  1): 
one  following  the  Herald  Canyon  east  of  Wrangell  Island 
and  another  through  the  Barrow  Canyon  between  the 
Alaskan  coast  and  Hanna  Shoals.  In  the  present  study, 
transports  through  both  pathways  are  selected  [among 
the  transports  through  other  sections  and  the  Bering 
Strait  transport  (BST)  itself]  as  TQs  used  for  optimizing 
the  HFR  positions. 

Prior  information  of  the  local  currents  is  specified  by 
assimilating  the  available  climatological  data  into  a  nu¬ 
merical  model.  Analysis  of  the  dynamical  balance  in  the 
assimilated  model  solution  made  by  Panteleev  et  al.  (2010) 
indicates  that  circulation  variability  in  the  BS  region  is 
governed  by  the  wind-forced  dynamics  of  a  topographically 
controlled  and  geostrophically  balanced  barotropic  flow. 

This  observation  justifies  monitoring  the  BST  using 
HFR  observations  of  surface  currents  and  thus  motivates 


Fig.  1.  Bathymetry  of  the  BS  region  (pink  contours,  m)  and  the 
major  pathways  of  Pacific  waters  (thick  blue  arrows).  Dashed  arrows 
show  alternating  coastal  current  north  of  the  Chukchi  Peninsula. 
Locations  of  the  tested  HFR  sites  (circles),  sectors  monitored  by  the 
HFRs,  and  their  coverage  by  observation  points  (for  Shishmareff 
location)  are  shown  in  red.  Straight  black  lines  denote  MHS  trans¬ 
port  sections  used  for  optimization  of  the  HFR  sites.  Locations  of 
the  permanent  NOAA  moorings  are  shown  by  yellow  diamonds. 
Boundary  of  the  numerical  model  domain  is  given  by  the  gray 
rectangle. 

our  optimization  of  HFR  locations  through  dynamically 
constrained  data  assimilation  combined  with  sensitivity 
analysis. 

Such  an  approach  for  configuring  long-term  HFR  ob¬ 
servations  may  benefit  the  overall  accuracy  of  measuring 
the  exchange  between  the  Arctic  and  Pacific  Oceans 
under  the  constraints  imposed  by  the  limited  number  of 
logistically  feasible  monitoring  sites  in  the  region.  True 
optimization  of  an  HFR  network  involves  many  factors 
and  guidelines  that  are  difficult  to  express  in  mathemat¬ 
ical  terms.  In  the  present  study  we  limit  ourselves  by  two 
principles:  1)  the  best  (optimal)  configuration  maximizes 
the  error  reduction  of  a  given  set  of  TQs  and  2)  logisti¬ 
cally  permissible  HFR  sites  are  limited  to  the  inhabited 
settlements  on  the  Alaskan  coast.  Of  course,  there  are 
many  more  guidelines  and  restrictions  imposed  by  envi¬ 
ronmental,  logistic,  and  financial  conditions,  and  it  is 
difficult  to  permanently  define  the  set  of  TQs  that  will 
satisfy  all  the  needs  of  the  scientific,  business,  and  local 
communities.  Nevertheless,  we  assume  that  the  proposed 
technique  is  worth  considering  as  an  element  in  the  de¬ 
sign  of  observational  networks  in  the  Arctic  and  in  the 
Bering  Strait  in  particular. 

Historically,  the  BST  has  been  estimated  from  moored 
observations  conducted  as  early  as  the  1930s  (Coachman 
and  Aagaard  1966,  1988;  Coachman  et  al.  1975),  but  its 
continuous  monitoring  began  in  1990  with  the  deployment 
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of  a  mooring  array  as  a  part  of  a  joint  program  between 
American  and  Russian  scientists  (Roach  et  al.  1995).  Since 
that  time,  at  least  three  moorings  have  been  maintained  in 
the  BS  region  (Fig.  1),  which  provide  a  20-yr  series  of  the 
BST  and  associated  MHS  fluxes  (Woodgate  et  al.  2005, 
2012).  In  that  context,  it  is  instructive  to  compare  the  ef¬ 
ficiency  of  monitoring  the  BST  by  a  limited  number  of 
moorings  and  by  a  pair  of  optimally  positioned  HFRs. 

Mathematically,  the  proposed  optimization  technique 
is  applicable  under  the  assumption  that  perturbations  of 
the  optimal  state  introduced  by  new  observations  are 
either  small  or  that  the  state  is  governed  by  quasi-linear 
dynamics.  For  the  large-scale  circulation  in  the  BS  re¬ 
gion,  both  assumptions  appear  to  be  acceptable  for  two 
reasons:  First,  as  it  was  shown  in  Panteleev  et  al.  (2010), 
the  nonlinear  terms  in  the  governing  equations  are  small 
compared  to  the  wind  forcing  and  geostrophy;  second, 
the  HFR  observation  errors  (which  include  model  rep¬ 
resentation  errors)  are  larger  than  the  respective  surface 
velocity  errors  determined  from  the  error  analysis  of  the 
optimal  state,  which  is  constrained  by  massive  in  situ  and 
remotely  sensed  observations. 

In  this  study,  the  data-optimized  climatological  evolu¬ 
tion  of  the  ocean  state  in  the  BS  region  during  the  ice-free 
period  (July-September)  is  used  as  a  reference  for  opti¬ 
mization  of  HFR  locations.  The  optimization  criterion  is 
the  overall  error  reduction  of  the  MHS  transport  estimates 
through  a  number  of  sections  in  the  southern  Chukchi  and 
northern  Bering  Seas,  including  the  Bering  Strait.  In  the 
next  section,  we  briefly  describe  the  basic  principles  of 
optimizing  the  observations,  the  numerical  model,  and  the 
methodology  of  the  numerical  experiments,  including 
a  description  of  the  observing  system  simulation  experi¬ 
ments  (OSSEs)  used  to  verify  the  results.  Section  3  con¬ 
tains  the  results  of  HFR  optimization  together  with 
recommendations  for  positioning  of  the  HFRs  on  the  U.S. 
coastline  and  comparison  of  the  monitoring  efficiencies  of 
the  BST  by  HFRs  and  by  permanent  moorings.  A  sum¬ 
mary  and  a  discussion  in  section  4  complete  the  paper. 

2.  Methodology 

a.  Optimizing  observations  under  the  dynamical 
constraints 

Consider  an  ocean  general  circulation  model  (OGCM) 
as  a  nonlinear  mapping  M  of  a  vector  of  control  param¬ 
eters  c  to  the  ocean  state  x: 

x  =  A/'(c).  (1) 

Here,  the  Af-dimensional  ocean  state  vector  x  e  1ZN  is 
composed  of  all  the  OGCM  fields’  values  discretized  on 
the  four-dimensional  space-time  grid,  while  the  control 


vector  c  contains  poorly  known  input  parameters  of  the 
model,  which  may  include  gridded  fields  of  the  initial 
and  open  boundary  conditions,  and  the  forcing  fields. 

Let  x0  provide  a  minimum  to  some  quadratic  cost  func¬ 
tion  /(x)  and  8x  be  small  deviations  from  xa.  Then  a  varia¬ 
tion  8J  of  J  from  its  value  at  the  minimum  is 

SJ  =  ^SxtB_1Sx  =  ^SctMtB^1MSc,  (2) 

where  B-1  is  a  nonnegative  definite  matrix  and  the 
superscript T  stands  for  transposition.  The  symbol  M  here 
is  the  linearization  of  OGCM  operator  J\f  in  the  vicinity 
of  cQ ,  the  control  vector  mapped  by  M  to  x0,  so  that 

8x  =  MSc  (3) 

holds. 

The  cost  function  [Eq.  (2)]  can  be  interpreted  as  the 
argument  of  a  Gaussian  probability  density  function.  In 
this  sense,  x0  can  be  referred  to  as  the  “optimal”  (most 
probable)  state  and  B-1  is  the  inverse  error  covariance 
matrix. 

Assume  that  K  additional  observations  d  of  the  state 
vector  x  are  made  such  that  d  =  Dx  +  e,  where  e  is  the 
observational  noise  with  covariance  W  and  D  is  the  linear 
operator  mapping  of  the  state  to  observations.  Given 
these  new  data,  the  variation  of  the  cost  function  [Eq.  (2)] 
should  be  augmented  with  the  additional  term 

SJ  ->  SJ  +  i(DSx  +  Sd)TW_1(DSx  +  Sd) ,  (4) 

where  Sd  =  Dx0  -  d  and  Sx  is  the  perturbation  of  the 
optimal  state  caused  by  the  cost  function  update.  It 
follows  from  Eqs.  (2)  and  (4)  that  the  updated  inverse 
error  covariance  is 

B_1  =  B_1  +  DtW_1D.  (5) 

Using  the  Woodbury  identity,  Eq.  (5)  can  be  rewritten 
in  the  form 

B  =  B-BDtW_1DB,  (6) 

where 

W  =  W  +  DBDt.  (7) 

Since  both  terms  in  the  rhs  of  Eq.(6)  are  nonnegative  def¬ 
inite,  adding  new  observations  always  reduces  the  overall 
error  variance  described  by  the  trace  of  B. 

In  oceanographic  practice,  observational  arrays  are 
often  designed  for  monitoring  certain  quantities,  such  as 
MHS  transports  across  key  sections.  By  assuming  that 
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a  set  of  m  target  quantities  (TQs)  q  =  {qi,  . . .  ,  qm}  is 
linearly  related  to  the  ocean  state, 

q  =  Qx,  (8) 

and  taking  Eq.  (6)  into  the  account,  we  can  deduce  that 
the  total  error  variance  of  q, 

B9  =  QBQt,  (9) 


is  reduced  by 

trSB^  =  tr(B?  -  B?)  =  tr(QBDTVir1DBQT)  =  tr(STS). 

(10) 


Here, 

s  =  w  1/2DBQt  (11) 

is  the  sensitivity  matrix  (e.g.,  Kohl  and  Stammer 
2004).  By  electing  the  trace  of  the  covariance  matrix  as 
its  norm,  we  employ  the  so-called  A-optimality  ap¬ 
proach  (e.g.,  Alvarez  and  Mourre  2012)  commonly 
used  in  geosciences  for  optimal  design  of  observa¬ 
tional  networks. 

If  D  is  the  N  X  N  identity  matrix  (i.e.,  full,  direct  state 
observation  at  all  considered  times),  then  the  columns 
of  S  span  the  subspace  of  7 ZN,  whose  elements  affect 
the  optimal  values  of  q.  In  particular,  if  m  =  1,  then 
there  is  a  single  TQ  (e.g.,  transport  through  a  section), 
and  the  sensitivity  matrix  is  given  by  an  A-dimensional 
vector  s.  Since  its  square,  sTs,  gives  the  reduction  of  the 
error  variance  [Eq.  (10)],  plotting  the  4D  distribution 
of  the  absolute  values  of  the  components  of  s  shows  the 
relative  contribution  of  the  observations  at  various  grid 
points  to  the  reduction  of  the  a  posteriori  error  vari¬ 
ance  of  the  TQ.  Therefore,  the  space-time  locations  of 
in  situ  observations  with  the  largest  impact  on  error 
variance  reduction  may  be  identified  by  inspecting  the 
columns  of  S.  To  simplify  the  visualization  of  the  4D 
sensitivities,  their  absolute  values  can  be  averaged  over 
certain  dimensions  (e.g.,  over  time)  to  produce  “sen¬ 
sitivity  maps”  that  show  the  relative  contributions  of 
various  observation  points  in  space  to  the  overall  error 
reduction  (see  section  3a). 

In  the  considered  case  of  HFR  observations,  a  single 
radar  measures  projections  of  the  surface  currents  on  its 
beams  at  hundreds  of  locations  (gray  dots  in  Fig.  1),  so 
that  the  operator  D  projects  gridded  velocity  vectors  of 
the  model  onto  the  beam  directions  at  all  K  observation 
points  of  the  radar.  A  two-radar  column  of  the  sensitivity 
matrix  is  therefore  K\  +  K2  elements  long,  where  Ki  and 
K2  are  the  numbers  of  observation  points  monitored  by 
the  individual  radars. 


From  the  computational  point  of  view,  the  most  ex¬ 
pensive  part  is  obtaining  the  columns  of  the  sensitivity 
matrix,  mostly  because  of  the  huge  size  of  B.  Furthermore, 
computing  the  action  of  B  on  a  column  q  of  QT  in  Eq.  (11) 
is  not  straightforward  because  the  cost  function  [Eq.  (2)]  is 
formulated  in  terms  of  the  sparse  rank-deficient  matrix 
B-1.  This  action  may  be  calculated  by  inverting  B-1  in  the 
subspace  spanned  by  control  variable  variations  Sc.  Using 
Bc  to  denote  this  inverse  (i.e.,  the  covariance  ScScT)  and 
taking  Eq.  (3)  into  the  account,  the  relationship 

B  =  SxSx r  =  MSc(MSc)t  =  MBcMt  (12) 

is  established,  where  overlines  [( )]  denote  the  ensemble 
average.  Term  Bc  has  a  significantly  smaller  dimension 
than  B  and  its  action  on  a  control  vector  can  be  com¬ 
puted  using  the  4DVar  code.  Substitution  of  Eq.  (12) 
into  Eq.  (11)  yields  the  basic  formula  for  computing  the 
sensitivity  matrix: 

S  =  W“1/2DMBcMtQt.  (13) 

To  speed  up  the  computations,  Bc  was  approximated 
by  the  low-rank  error  covariance  matrix  retrieved  from 
the  data  assimilative  run  of  the  Arctic  Cap  Hybrid  Co¬ 
ordinate  Ocean  Model  (AC/HYCOM)  (section  2c). 

The  matrix  M  is  implicitly  given  by  the  OGCM  code 
linearized  in  the  vicinity  of  the  control  vector  associated 
with  the  model  optimal  state,  both  described  in  the  next 
section.  The  code  for  computing  the  action  of  MT  on 
a  vector  was  available  in  the  form  of  the  adjoint  model. 
The  TQs  were  given  by  the  set  of  the  mean  MHS  fluxes 
through  eight  sections  shown  in  Fig.  1.  Mathematically, 
these  fluxes  are  described  by  the  integrals  over  the  target 
sections  at  with  /  =  1 . m  as  follows: 

(14) 


v(6  —  0q)  daj ; 

S-S0)daj,  (15) 

where  angular  brackets  ((•))  denote  time  average,  v  is  the 
velocity  component  perpendicular  to  the  ith  section,  6  is 
the  seawater  temperature,  S  is  salinity,  and  Cp  is  the 
specific  heat,  where  60  =  —  1.9°C  and  S0  =  34.8  psu  are  the 
background  values  used  for  computing  HS  fluxes  through 
the  Bering  Strait  (e.g.,  Woodgate  et  al.  2012). 

The  rows  q,  of  the  TQ  matrix  Q  in  Eq.  (8)  were  com¬ 
puted  by  linearizing  Eq.  (15)  in  the  vicinity  of  the  refer¬ 
ence  state  and  then  differentiating  the  finite-difference 
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analogs  of  Eqs.  (14)  and  (15)  with  respect  to  all  the 
gridpoint  values  of  v,  0,  and  S. 

As  a  result,  a  column  of  S  was  computed  in  five  steps: 

1)  Compute  q'  =  MTq  by  running  the  adjoint  model 
forced  by  q. 

2)  Compute  c  =  Bcq'  by  convolving  the  error  covari¬ 
ance  Bc  between  the  control  variables  of  the  model 
with  q'. 

3)  Compute  x  =  Me  by  running  the  tangent  linear  model 
forced  by  c. 

4)  Project  x  onto  the  observational  points  of  a  given  pair 
of  HFRs  by  computing  the  product  x  =  Dx. 

5)  Rescale  the  result  using  multiplication  by  the  inverse 
of  the  HFR  measurement  error:  s  =  W_1/2x. 

In  the  following  sections,  we  give  more  detail  on  the 
model  M  and  estimations  of  the  reference  state  and  error 
covariance  matrices  Bc  and  W. 

b.  Background  state  in  the  BS  region 

The  reference  state  was  found  by  assimilating  the 
monthly  averaged  oceanographic  (SST,  temperature, 
salinity  velocity,  and  ice  observations)  and  atmospheric 
(NCEP  reanalysis)  data  into  the  semi-implicit  ocean 
model  (SIOM).  The  monthly  averaged  temperature- 
salinity  (T-S)  data  and  corresponding  variances  were 
obtained  by  averaging  over  the  model  grid  cells  all  of  the 
available  17  960  temperature  and  13  380  salinity  profiles 
for  July-September.  The  number  of  T-S  profiles  in  some 
model  grid  boxes  was  insufficient  to  derive  reliable  es¬ 
timates  of  data  variance.  In  these  cases,  variability  was 
represented  by  horizontal  interpolations  of  statistics 
characterizing  data  variability.  The  obtained  climato¬ 
logical  data  and  the  estimates  of  the  corresponding 
standard  deviations  were  used  in  the  data  assimilation 
procedure.  The  estimates  of  the  monthly  transports 
through  the  Bering  Strait  were  taken  from  Woodgate 
et  al.  (2005). 

Ice  concentration  from  the  NCEP  reanalysis  was  used 
to  set  the  model  SST  to  the  freezing  value  if  the  ice 
concentration  exceeded  the  threshold  value  of  60%. 
NCEP-NCAR  wind  stress  and  surface  heat-salt  flux 
climatology  (http://www.cdc.noaa.gov/cdc/data.ncep. 
reanalysis.derived.html)  are  too  smooth  for  8-km  res¬ 
olution.  To  allow  for  the  adjustment  of  the  spatial  details 
in  the  model  forcing,  we  used  wind  stress  and  heat-salt 
flux  data  with  relatively  high  error  variance,  equal  to 
40%  of  their  spatial  and  temporal  variability  in  the 
northern  BS  and  Chukchi  Sea.  This  error  level  is  con¬ 
sistent  with  the  errors  in  the  NCEP-NCAR  forcing  no¬ 
ticed  by  Ladd  and  Bond  (2002). 

The  assimilation  procedure  is  similar  to  the  one  used 
by  Panteleev  et  al.  (2007)  to  study  the  climatological 


circulation  in  the  Kara  Sea,  and  by  Luchin  and  Panteleev 
(2014)  to  study  climatological  states  in  the  Chukchi 
Sea,  with  the  only  difference  being  that  the  assimi¬ 
lated  monthly  climatologies  (July-September)  were 
derived  for  the  BS  region  (Fig.  1)  where  the  model  was 
configured. 

The  SIOM  is  a  modification  of  the  C-grid,  z-coordinate 
OGCM  developed  in  the  Laboratoire  d’Oceanographie 
Dynamique  et  de  Climatologie  (Madec  et  al.  1998).  The 
model  was  specifically  designed  for  the  implementa¬ 
tion  of  4DVar  methods  into  regional  simulations  con¬ 
trolled  by  fluxes  at  the  open  boundaries  and  sea  surface. 
The  OGCM  is  semi-implicit  with  respect  to  both  baro- 
tropic  and  baroclinic  modes,  permitting  simulations 
with  relatively  large  time  steps  of  approximately  0.1 
day  (Nechaev  et  al.  2005;  Panteleev  et  al.  2006a, b).  The 
tangent  linear  model  was  obtained  by  direct  differenti¬ 
ation  of  the  forward  model  code.  The  adjoint  code  of 
the  model  was  built  analytically  by  transposition  of  the 
tangent  linear  model  operator  linearized  in  the  vicinity 
of  the  given  solution  of  the  forward  model.  Further 
details  of  the  SIOM’s  numerical  scheme  can  be  found  in 
Nechaev  et  al.  (2005). 

In  the  present  study,  the  computational  grid  has 
a  horizontal  resolution  of  8  km,  which  is  fine  to  resolve 
seasonal  dynamics  from  a  climatological  perspective  and 
coarse  enough  to  suppress  instabilities  of  the  adjoint 
model  at  scales  comparable  with  the  local  deformation 
radius  (7  km).  This  permits  use  of  a  2-month-long  as¬ 
similation  interval  (from  16  July  to  15  September)  for 
integrating  the  model  and  its  adjoint. 

In  the  vertical  dimension,  the  model  has  13  unequally 
spaced  layers  ranging  in  thickness  from  2.5  m  at  the  surface 
to  10  m  near  the  bottom.  The  maximum  depth  of  the  re¬ 
gion  is  86  m  (Fig.  1).  The  model  is  controlled  by  the  initial 
fields  of  temperature,  salinity,  velocity,  and  surface  ele¬ 
vation  (on  16  July),  and  by  the  three  monthly-mean  values 
of  the  same  fields  at  the  open  boundaries.  The  number  of 
gridded  values  of  the  control  variables  at  the  start  of  in¬ 
tegration  was  144397.  With  the  number  of  the  gridded 
values  of  temperature,  salinity,  and  velocity  at  the  open 
boundaries  equal  to  13389,  the  total  number  of  control 
variables  (dimension  n  of  the  control  vector  c)  is  157  786. 
The  number  of  assimilated  observations  exceeded  the 
number  of  control  variables  several  times,  rendering  the 
assimilation  problem  relatively  well  determined. 

Evolution  of  the  reference  state  was  obtained  by 
finding  the  values  of  control  variables  that  minimized 
the  distance  between  the  model  fields  and  observations. 
The  latter  were  specified  once  a  day  (every  tenth  time 
step)  by  linear  interpolation  of  the  three  climatologies  in 
time.  With  the  total  duration  of  62  days  (from  the  be¬ 
ginning  of  16  July  to  the  end  of  15  September),  there 
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were  63  time  steps  when  observations  were  available. 
Fields  at  the  model  boundaries  were  specified  in  a  simi¬ 
lar  manner  by  linear  interpolation  in  time  of  the  control 
values  on  day  1  (beginning  of  16  July),  day  31  (16  Au¬ 
gust),  and  day  62  (end  of  15  September).  Before  the 
assimilation,  these  values  were  taken  from  the  clima¬ 
tologies  and  then  varied  together  with  the  initial  fields 
(on  16  July)  in  the  process  of  optimization. 

As  an  example,  evolution  of  the  optimized  sea  surface 
temperature  (SST)  and  velocity  fields  is  shown  in  Fig.  2. 
The  rms  difference  of  the  optimized  temperature  fields 
from  the  July  and  September  climatologies  (0.5°-0.9°C) 
lies  well  within  the  error  bars.  In  addition,  the  optimized 
fields  are  constrained  by  the  dynamics  of  the  SIOM. 
Similar  correspondence  was  found  for  the  monthly-mean 
salinity  fields.  The  evolution  of  the  optimized  Bering 
Strait  transport  (Fig.  3)  shows  the  well-documented  trend 
of  the  BST  seasonal  cycle  (Woodgate  et  al.  2005,  2010, 
2012),  and  it  also  fits  the  observational  error  bars  (0.3  Sv) 
of  Woodgate  et  al.  (2005),  shown  by  the  vertical  sides  of 
the  gray  rectangles  in  Fig.  3.  Optimized  monthly-mean 
transport  of  freshwater  and  heat  through  the  Bering 
Strait  also  demonstrate  a  similar  consistence  with  the 
estimates  of  Woodgate  et  al.  (2005),  primarily  because 
the  model  is  constrained  by  the  massive  hydrological 
dataset  (17  960  temperature  and  13  380  salinity  profiles) 
employed  in  the  assimilation  process. 

Optimized  transports  through  the  other  key  sections 
(shown  by  straight  black  lines  in  Fig.  1)  have  also  been 
checked  against  (rather  scanty)  estimates  available  in 
literature:  The  optimal  summer  transport  of  the  Anadyr 
Current  between  the  Chukchi  Peninsula  and  St.  Lawrence 
Island  is  0.8  ±  0.2  Sv  (where  lSv  =  106m3s_1)  and  is  in 
good  agreement  with  the  estimate  of  Coachman  (1993), 
who  obtained  the  value  of  0.7  ±  0.2  Sv  from  1985  to  1988 
current  meter  data. 

The  optimized  flow  pattern  in  the  northern  part  of 
the  model  domain  (Fig.  2c)  agrees  well  with  known 
schemes  of  the  regional  circulation  (e.g.,  Winsor  and 
Chapman  2004;  Clement  et  al.  2005),  revealing  the 
central  northward  flow  and  flow  through  the  Gulf  of 
Kotzebue  (solid  blue  arrows  in  Fig.  1).  The  ratio  be¬ 
tween  the  optimized  transports  of  these  branches 
(0.34  Sv/0.47  Sv)  is  in  a  good  agreement  with  the  ratio 
of  the  annual  mean  transports  through  the  Herald 
Channel  and  the  section  between  the  Herald  Shoal  and 
Cape  Lisburne  obtained  by  Panteleev  et  al.  (2010).  The 
optimized  circulation  also  reveals  the  well-known  weak 
(—0.1  Sv)  East  Siberian  Current  (ESC)  along  the  Siberian 
Coast  (Weingartner  et  al.  1999). 

Based  on  these  comparisons  against  independent 
observations,  we  assume  that  the  optimized  solution 
realistically  simulates  the  major  features  of  the  local 


Fig.  2.  SST  (°C)  of  the  reference  state  in  (a)  July  and  (b)  September; 
and  (c)  the  respective  seasonal  mean  SSH  (cm)  and  surface  velocity 
fields.  Every  second  arrow  is  shown. 
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Fig.  3.  Evolution  of  the  AC/HYCOM  BST  averaged  over  two 
seasons — 2011  and  2012  (thick  solid  line).  Dashed  lines  are  the 
monthly-mean  values.  Evolution  of  the  monthly-mean  SIOM  trans¬ 
ports  is  given  by  thin  solid  lines.  Gray  rectangles  show  the  respective 
climatological  estimates  of  Woodgate  et  al.  (2005). 

circulation  in  July-September  and  thus  gives  us 
grounds  to  use  it  as  a  background  state  x0  in  the  search 
for  an  optimal  configuration  of  HFR  sites. 

c.  Estimation  of  the  error  covariances 

To  assess  the  error  covariance  of  the  reference  state  [Bc 
in  Eq.  (13)],  we  analyzed  a  multiyear  run  of  an  OGCM 
sequentially  constrained  by  the  data.  For  this  purpose  we 
used  daily  analyses  together  with  the  12-  and  24-h  forecasts 
of  the  V120  AC/HYCOM  (Posey  et  al.  2010;  Allard  et  al. 
2012)  runs  in  July-September  of  2011  and  2012.  The  error 
statistics  was  collected  from  the  ensemble  of  differences 
between  the  AC/HYCOM  analyses  and  the  forecasts 
projected  onto  the  SIOM  grid.  As  a  result,  the  ensemble 
had  k  =  90  days  X  2  =  180  daily  sampled  members.  Since 
an  ensemble  member  provides  a  sample  of  the  error  fields 
of  the  control  variables,  the  sample  error  covariance  in  the 
control  space  Bc  can  be  computed  by 

Br-AiCCT'  <16) 

where  C  is  the  k  X  n  matrix  of  the  column-wise  listed 
ensemble  members  with  the  ensemble  means  removed. 
As  it  is  seen  from  Eq.  (16),  the  nXn  matrix  Bc  is  strongly 
degenerate,  as  it  has  no  more  than  k  -  1  <C  n  non-0  ei¬ 
genvalues.  However,  its  product  with  a  control  vector 
requires  only  2 kn  <  n 2  multiplications.  To  increase  the 
computational  efficiency,  we  performed  a  singular  value 
decomposition  of  Bc,  which  also  gives  insight  into  the 
spatial  structure  of  the  dominant  error  modes. 

The  spectrum  of  Bc  computed  with  respect  to  the 
energy  norm  shows  that  the  first  eight  eigenvectors  ac¬ 
count  for  90%  of  the  error  variability.  In  the  actual 
computations  involving  multiplication  by  Bc  (e.g.,  step  2, 


Fig.  4.  Sea  surface  height  (cm)  and  velocity  at  z  =  5  m  of  the  first 
eigenvector  of  the  error  covariance  between  the  control  variables 
corresponding  to  the  initial  conditions. 


section  2a),  we  used  10  eigenvectors  accounting  for  93% 
of  the  total  error  estimated  from  the  ensemble. 

It  is  also  instructive  to  analyze  the  leading  eigenvectors 
of  Bc.  Figures  4  and  5  show  parts  of  the  leading  eigen¬ 
vector  that  correspond  to  the  initial  conditions  in  sea 
surface  height  and  velocity  (Fig.  4),  and  the  August  mean 
values  of  the  meridional  velocity  and  temperature  at  the 
northern  open  boundary  (Fig.  5).  Figure  4  indicates  that 
the  leading  error  mode  is  close  to  geostrophic  balance 
with  somewhat  larger  uncertainties  in  the  velocities  along 
the  U.S.  side  of  the  Bering  Strait.  Figure  5  demonstrates 
that  errors  in  the  Alaska  Coastal  Current  (west-most  part 
of  the  section)  and  the  Herald  Canyon  flow  (172°-173°W, 
see  Fig.  1)  tend  to  be  negatively  correlated.  Temperature 
errors  seem  to  be  dynamically  insignificant;  supporting 
the  notion  that  circulation  in  the  region  is  predominantly 
barotropic  in  nature. 

With  the  error  covariance  model  [Eq.  (16)]  at  hand, 
estimation  of  the  second  term  in  W  =  W  +  DBDT  [Eq. 
(7)]  also  becomes  computationally  feasible.  We  found 
that  the  rms  error  variances  of  the  radial  velocities  ex¬ 
tracted  from  the  reference  state  varied  within  2-4  cms"1 
and  therefore  can  be  neglected  compared  to  the  HFR 
observational  error  variances  W,  which  have  a  typical 
magnitude  of  5-15 cms-1.  Taking  this  into  account,  the 
final  approximation  for  the  sensitivity  matrix  (13)  has 
the  form 

s  =  c,;iv1dmbmtqt,  (17) 

where  V  is  the  (diagonal)  matrix  of  the  HFR  rms  error 
variances  and  Cr  is  the  square  root  of  the  HFR  error 
correlation  matrix. 
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Fig.  5.  Meridional  velocity  (black  contours,  cms-1)  and  temperature  (white  contours,  0.01°C) 
components  of  the  first  eigenvector  of  the  error  covariance  between  the  control  variables  at  the 
northern  open  boundary  of  the  computational  domain. 


The  error  covariance  model  [Eq.  (16)]  also  simplifies 
the  estimation  of  the  background  error  covariance 
matrix  Bq  [Eq.  (9)]  between  the  TQs.  The  magnitude  of 
its  trace  provides  a  reference  for  assessing  the  relative 
rms  error  reduction  in  a  specified  set  of  target 
quantities  Q: 


ArSB  \1/2 

’  tr(STS)  " 

ltrBJ 

tr(QBQT) 

In  this  study,  the  three  sets — QM,  QH ,  and  — of 
target  quantities  are  specified  by  the  MHS  transports 
through  the  sections  shown  in  Fig.  1. 

d.  HFR  observations 

Since  the  TQs  used  for  optimization  of  HFR  posi¬ 
tions  are  characterized  by  a  scale  typically  on  the  order 
of  tens  of  kilometers,  the  HFR  observational  operators 
were  chosen  to  mimic  the  long-range  5-MHz  Coastal 
Ocean  Dynamics  Application  Radar  (CODAR)  sys¬ 
tem  (Barrick  et  al.  1977).  This  type  of  HFR  employs 
two  crossed-loop  antennas  around  a  whip  (vertical 
monopole)  as  a  receiving  antenna  and  two  whips  as 
a  transmitting  antenna,  and  is  capable  of  observing  sea 
surface  velocities  up  to  250  km  offshore  at  radial,  an¬ 
gular,  and  temporal  resolutions  of  6  km,  5°,  and  20  min, 
respectively  (Fig.  1). 

Recent  estimates  of  the  CODAR  radial  velocity  er¬ 
rors  (e.g.,  Lipa  2003;  Emery  et  al.  2004;  Yoshikawa  et  al. 
2006;  Laws  et  al.  2010)  indicate  their  variation  from  3  to 
4 cms-1  at  small  ranges  to  15-20 cms-1  at  200-250 km. 
Based  on  this,  we  assumed  a  linear  increase  of  the  rms 
error  variances  with  distance  from  the  HFR  base  and 
adopted  the  following  analytical  model  for  the  diagonal 
elements  of  the  observation  error  variance  matrix  V  in 
Eq.  (17): 


Vi  =  v0(l  +  r.fr0),  (19) 

where  rt  is  the  distance  of  the  ith  observation  point  from 
the  parent  radar,  v0  =  3cms_1,  and  r0  =  50  km. 

The  direction-finding  method  (Lipa  and  Barrick  1983) 
employed  by  the  CODAR  system  results  in  certain  er¬ 
rors  (3°-10°)  in  the  bearing  angle  of  the  received  signal. 
This  property  induces  correlations  between  the  radial 
velocity  errors  of  the  adjacent  beams.  To  account  for 
this,  the  entries  of  Cr  were  defined  by 

Cy  =  exp( — d/r1 )  cos (ad) .  (20) 

Here  d  is  the  distance  between  two  observation  points 
(not  necessarily  monitored  by  the  same  radar),  ad  is  the 
angle  between  the  radial  velocities  observed  in  these 
points,  and  u  =  10  km. 

In  the  experiments  below,  we  consider  only  pairs  of 
HFRs  because  single  coastal  radars  are  rarely  deployed: 
the  major  purpose  of  such  deployments  is  monitoring  and 
prediction  of  the  surface  transport,  which  has  many  prac¬ 
tical  applications  (search  and  rescue,  pollution,  etc.)  and 
cannot  be  effectively  carried  out  by  observing  only  one 
velocity  component.  Considering  more  than  two  radars 
would  require  more  voluminous  computations,  which,  in 
principle,  could  be  carried  out,  but  lie  beyond  the  scope  of 
the  present  study,  whose  major  objective  is  to  demonstrate 
the  performance  of  the  sensitivity  technique  in  application 
to  the  HFR  deployment. 

3.  Results 

a.  Impact  of  the  surface  velocity  observations 

Figure  6  demonstrates  the  time-averaged  sensitivity 
(section  2a)  of  the  mean  Bering  Strait  mass  transport  to 
observations  of  surface  velocities  at  the  grid  points  of  the 
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Fig.  6.  Time-averaged  sensitivity  map  of  the  mean  BST  to  sur¬ 
face  velocity  observations.  Sensitivity  values  are  normalized  by 
their  maximum. 

computational  domain.  The  distribution  was  computed  by 
setting  D  in  Eq.  (11)  to  have  unit  diagonal  values  for  sur¬ 
face  velocity  components  (i.e.,  surface  velocity  field  ob¬ 
served  in  all  grid  points  at  all  times)  and  zero  otherwise. 
The  map  provides  a  qualitative  insight  for  optimally  po¬ 
sitioning  a  pair  of  HFRs  among  the  seven  locations 
available  within  the  existing  U.S.  settlements  in  the  region 
(Fig.  1).  Larger  sensitivity  values  on  the  map  correspond  to 
regions  where  surface  velocity  observations  provide 
a  larger  reduction  of  the  a  posteriori  error  variance  of  the 
BST.  Ideally,  a  pair  of  radars  should  be  positioned  to  cover 
the  maximum  possible  area  of  the  high-sensitivity  regions 
in  Fig.  6.  These  regions  tend  to  align  along  the  major 
currents  contributing  to  the  BST. 

Examination  of  the  maps  in  Figs.  1  and  3  shows  that 
HFR  placement  at  Cape  Prince  of  Wales  (CPW;  point  4) 
and  at  Little  Diomede  Island  (DI;  point  5)  are  the  best  to 
cover  the  hot  spots  in  the  vicinity  of  the  Bering  Strait, 
whereas  HFRs  at  Savonga  (point  1)  and  Gambell  (point 
3)  are  able  to  cover  the  high-sensitivity  patch  associated 
with  the  Anadyr  Current,  around  64.2°N  172.3°W.  HFR 
deployment  at  Shishmareff  (point  6)  or  Point  Hope 
(point  7)  appears  to  be  less  efficient  because  radars  at 
these  locations  are  able  to  accurately  monitor  only  the 
low-sensitivity  areas  (Fig.  6). 

From  the  “common  sense”  point  of  view,  one  may 
argue  that  placement  of  radars  close  to  the  BS  (at  DI  and 
CPW)  should  be  the  best  because  they  measure  the  BS 
flow  directly  at  the  ranges  of  30-60  km  (i.e.,  with  the 
highest  accuracy).  However,  one  should  keep  in  mind 
that  HFRs  measure  surface  currents  with  large  wave-  and 
wind-driven  components  that  are  less  correlated  with  the 


long-term  BST  variations  than  the  geostrophic  compo¬ 
nent.  In  that  respect  simultaneous  monitoring  of  the  BST 
and  the  Anadyr  Current  south  of  the  Chukchi  Peninsula 
may  provide  more  information  about  the  mean  BST  by 
way  of  the  geostrophically  induced  horizontal  correla¬ 
tions  between  the  surface  currents  at  64°-64.5°N,  170°- 
174°W  and  the  BST. 

Furthermore,  the  problem  of  optimal  positioning  be¬ 
comes  much  more  difficult  if  we  consider  a  large  set  of 
the  TQs,  such  as  MHS  transports  through  multiple  sec¬ 
tions  shown  in  Fig.  1.  In  such  a  complex  situation,  the 
optimization  approach  described  in  section  2  may  pro¬ 
vide  useful  guidance  for  configuring  HFR  networks. 

b.  Optimal  HFR  locations 

Using  the  algorithm  described  in  section  2a,  a  series  of 
numerical  experiments  with  pairs  of  HFRs  virtually 
deployed  in  the  seven  settlements  shown  in  Fig.  1  were 
conducted.  All  HFRs  had  identical  characteristics:  the 
aperture  of  150°  and  the  range  of  200  km  with  the  middle 
beam  oriented  perpendicular  to  the  coastline.  The  only 
exceptions  were  the  CPW  (point  4)  and  DI  (point  5) 
locations,  where  the  northward  and  southward  orienta¬ 
tions  of  the  middle  beam  were  considered.  The  obser¬ 
vation  errors  of  the  radars  were  assumed  to  be  described 
by  Eqs.  (19)  and  (20).  The  azimuthal  and  radial  resolu¬ 
tions  were  5°  and  6  km,  respectively. 

Figure  7  shows  the  values  of  the  relative  error 
reductions — and  ^ — for  36  different  combina¬ 
tions  of  the  HFR  pairs.  The  thick  gray  curve  for  the  mass 
transport  variance  reduction  shows  two  distinctive  peaks. 
The  first  one  involves  a  pair  of  southward-looking  HFRs 
positioned  at  CPW  (point  4)  and  DI  (point  5),  capable  of 
reducing  errors  in  the  mass  transports  by  6.3%.  The 
second  peak  (6. 0-5. 9%  reduction)  corresponds  to  the 
HFR  deployments  in  Sinuk  (point  2)  and  either  DI  or 
CPW.  This  configuration  provides  more  accurate  moni¬ 
toring  of  the  U.S.  part  of  the  BST,  which  contributes 
70%-80%  to  the  total  transport  (Coachman  et  al.  1975; 
Coachman  and  Aagaard  1988).  Overall,  optimal  HFR 
deployment  provides  1.3-3  times  better  reduction  of  the 
transport  errors  compared  to  other  configurations. 

It  is  noteworthy  that  pairs  involving  the  key  northward¬ 
looking  DI  and  CPW  radars  demonstrate  somewhat 
smaller  error  reduction  compared  to  southward-looking 
configurations.  The  increased  efficiency  of  the  southward¬ 
facing  HFR  is  due  to  the  longer  residence  time  of  observed 
information  as  it  propagates  through  the  computational 
domain,  which  allows  longer  residence  of  the  observed 
information  in  the  computational  domain.  Similar  effects 
have  been  also  obtained  in  the  efficiency  analysis  of  the 
hydrophysical  surveys  (e.g.,  Panteleev  and  Semyonov 
1988;  Panteleev  et  al.  2004). 
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HFR  pair 


Fig.  7.  Relative  reduction  £;q  of  the  errors  in  estimation  of  the  MHS  transports  for  various 
pairs  of  HFRs.  Numbers  labeling  HFR  pairs  correspond  to  locations  in  Fig.  1.  Bars  over  the 
CPW  (point  4)  and  DI  (point  5)  locations  denote  northward-looking  antennas  at  those  sites. 


It  is  also  evident  that  compared  to  the  volume  transport, 
the  curves  for  the  heat  and  salt  (HS)  transports’  error  re¬ 
duction  demonstrate  higher  peaks  for  the  Sinuk-DI  and 
Sinuk-CPW  pairs  (6.3%  and  6.1%,  respectively).  This  can 
be  explained  by  the  fact  that  HS  transports  in  July- 
September  are  particularly  sensitive  to  variations  in  the 
temperature  and  salinity  fields  of  the  eastern  part  of  the 
domain,  where  large  riverine  inflow  of  fresh  and  warm 
water  from  Alaska  drive  strong  horizontal  gradients  (Fig.  2). 

Pairs  involving  Shishmareff  (point  6)  and  Point  Hope 
(point  7)  locations  demonstrate  lower  values  of  i;q  in 
correspondence  with  the  qualitative  analysis  of  the 
sensitivity  map  in  the  previous  section.  The  overall  error 
reduction  (computed  by  averaging  the  three  curves  in 
Fig.  7)  gives  a  minor  preference  to  the  DI-CPW  pair 
(6.3%),  but  its  advantage  over  the  Sinuk-CPW  config¬ 
uration  (6.1%)  may  be  insignificant  given  the  un¬ 
certainties  in  the  model  and  observational  errors. 

Interestingly,  placement  of  the  oppositely  directed 
HFRs  at  DI  and  CPW  (44  and  55,  respectively,  config¬ 
urations  in  Fig.  7)  also  results  in  a  fairly  good  reduction 
of  error  (j;q  =  5.8%)  despite  the  fact  that  HFR  coverage 
areas  do  not  intersect  for  these  configurations.  The  re¬ 
sult  could  be  explained  by  the  fact  that  such  configura¬ 
tions  cover  a  much  larger  area  with  strong  sensitivities 
(Fig.  6)  and  weakly  correlated  (independent)  observa¬ 
tions.  Such  a  choice  of  observation  points  is  capable  of 
providing  a  fairly  large  amount  of  information  on  the 
magnitude  and  redistribution  of  the  BST  because 
the  measured  radial  velocities  tend  to  be  aligned  in  the 
principal  direction  of  the  flow. 

Although  the  relative  error  reduction  varies  2-3  times 
within  the  set  of  HFR  pairs,  the  maximum  magnitude  er¬ 
ror  reduction  may  seem  to  be  rather  small  (2%-6%)  given 
the  large  amount  of  surface  velocity  data  (N  ~  3  X  105 
observations)  produced  by  a  pair  of  radars  in  3  months. 
However,  the  value  of  N  is  still  relatively  small  compared 
to  the  number  of  observations  (—1.5  X  106)  used  for 


generating  the  reference  state.  Additionally,  these  obser¬ 
vations  cover  all  the  components  of  the  sea  state  (3D 
temperature,  salinity,  and  velocity)  and  therefore  formally 
define  the  seasonal  mean  values  of  the  TQs  with  an  ac¬ 
curacy  that  is  high  enough  to  restrict  HFR  data  from 
providing  a  significant  reduction  of  a  posteriori  error. 

In  that  respect,  it  is  instructive  to  compare  the  above- 
mentioned  values  of  i;q  with  the  values  provided  by  three 
moorings  (Woodgate  et  al.  2012)  permanently  main¬ 
tained  in  the  region  for  the  last  5  yr  (Fig.  1).  The  corre¬ 
sponding  computation  is  easily  performed  by  changing 
HFR  observational  operators  to  operators  that  sample 
three  vertical  profiles  of  temperature,  salinity,  and  ve¬ 
locity  from  the  model  fields  at  the  mooring  locations. 
The  resulting  values  of  %q  appeared  to  be  significantly 
(6-8  times)  lower  than  those  obtained  for  a  pair  of 
HFRs.  This  means  that  a  limited  number  of  moorings 
provide  much  less  additional  information  about  the  TQs 
than  a  pair  of  HFRs.  In  other  words,  the  key  circulation 
features  in  the  BS  region  could  be  better  monitored  by 
a  pair  of  HFRs.  This  conclusion  may  seem  intuitive 
because  HFRs  deliver  much  better  spatial  coverage  at 
approximately  the  same  temporal  resolution.  On  the 
other  hand,  moored  observational  platforms  provide 
information  in  the  vertical,  including  temperature  and 
salinity  profiles,  and,  more  importantly,  are  capable  of 
operating  under  ice  conditions. 

Inspecting  Fig.  7,  one  may  also  think  that  graphs  re¬ 
flect  intuitive  knowledge  that  is  often  used  in  planning 
oceanographic  observations.  We  assume,  however,  that 
rigorous  sensitivity  and  error  analysis  reflected  in 
Figs.  5-7  may  serve  as  a  reasonable  quantitative  guid¬ 
ance  in  planning  HFR  observations  in  the  Bering  Strait 
region.  For  example,  despite  a  somewhat  smaller  error 
reduction,  oppositely  directed  HFR  pairs  at  DI  or  CPW 
could  be  considered  as  reasonable  alternatives  to  the 
best  configurations  (DI-CPW  and  Sinuk-CPW)  because 
of  the  considerable  savings  in  their  maintenance  costs. 
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To  verify  the  above-mentioned  results,  a  series  of 
OSSEs  was  conducted. 

c.  Validation  of  the  results 

The  OSSE  technique  is  widely  used  for  skill  assessment 
of  the  data  assimilation  systems.  The  underlying  idea  is  to 
simulate  artificial  “data”  using  some  reference  model 
solution  as  a  “true”  ocean  state,  contaminate  these  data 
with  noise  (mimicking  observational  and  model  errors) 
to  produce  a  “first  guess”  solution,  and  then  try  to  re¬ 
construct  the  “true  state”  from  these  data. 

In  the  present  study,  the  OSSEs  were  used  to  validate 
the  optimized  HFR  locations  obtained  from  the  sensi¬ 
tivity  analysis.  “HFR  data”  were  extracted  from  the 
reference  solution  described  in  section  2b.  The  first -guess 
solutions  started  from  the  geostrophically  balanced  state 
corresponding  to  the  July  climatology  of  temperature  and 
salinity. 

Next,  simulated  HFR  data  for  a  tested  pair  of  loca¬ 
tions  were  assimilated  into  the  SIOM  using  the  4DVar 
algorithm  with  the  same  controls  and  data  (excluding  BS 
transport  observations)  that  were  used  for  optimization 
of  the  HFR  locations. 

The  following  model  for  the  simulated  HFR  obser¬ 
vations  v*  was  adopted: 

t>*  =  t>  +  CrVe,  (21) 

where  v  is  the  TGdimensional  vector  of  the  true  pro¬ 
jections  of  surface  velocity  on  the  HFR  beams  for 
a  given  pair  of  radars,  and  e  is  the  noise  from  the  stan¬ 
dard  normal  distribution. 

After  assimilating  the  HFR  data,  transports  through 
the  sections  in  Fig.  1  were  computed  and  the  root-mean- 
square  differences  between  the  reconstructed  transports 
and  their  values  from  the  reference  (true)  solution  were 
estimated.  These  mean  differences  and  the  correspond¬ 
ing  difference  for  the  BS  transport  were  used  as  quanti¬ 
tative  criterion  to  assess  the  quality  of  the  tested  pairwise 
HFR  locations.  Table  1  shows  the  results  of  the  OSS 
experiments  for  several  selected  HFR  pairs. 

As  it  is  seen  from  the  table,  there  is  a  fairly  good 
agreement  between  the  OSSEs  and  the  results  of  the 
sensitivity  analysis  at  Fig.  7:  the  best  HFR  configurations 
are  for  the  DI-CPW,  Sinuk-CPW,  and  for  the  oppo¬ 
sitely  directed  antennas  on  Little  Diomede  Island.  The 
Di-Point  Hope  and  Savonga-Sinuk  pairs  demonstrate 
significantly  larger  discrepancies,  with  the  worst  result 
provided  by  the  Savonga-Gambel  pair. 

Although  the  rms  difference  between  true  and 
reconstructed  volume  transports  can  be  considered  as 
a  reasonable  quantification  of  error  in  the  general  circu¬ 
lation  pattern,  accuracy  of  the  BS  transport  may  be  of 


Table  1.  Mean  rms  difference  between  the  true  and  assimilated 
transports  through  the  transects  shown  in  Fig.  1,  and  absolute  dif¬ 
ference  between  the  true  and  assimilated  BST  values.  Optimiza¬ 
tion  was  done  using  the  simulated  data  from  HFR  pairs. 


HFR  pair 

Mean  error  (Sv) 

BST  error 

Diomede-CPW 

0.092 

0.084 

Sinuk-CPW 

0.108 

0.129 

Diomede-Diomede 

0.123 

0.177 

Savonga-Sinuk 

0.192 

0.285 

Savonga-Gambel 

0.246 

0.454 

Diomede-Point  Hope 

0.169 

0.239 

particular  interest  because  of  its  strong  correlation  with  ice 
and  thermal  conditions  in  the  Chukchi  Sea  (Woodgate 
et  al.  2010;  Luchin  and  Panteleev  2014).  The  second  col¬ 
umn  of  the  table  shows  the  difference  between  re¬ 
constructed  and  true  transports  through  the  Bering  Strait. 
Taking  into  account  that  the  July-September  mean  BST  is 
about  1  Sv,  we  can  say  that  both  DI-CPW  and  Dl-Sinuk 
give  rather  accurate  estimates  of  the  BST,  while  posi¬ 
tioning  HFRs  in  Savonga-Sinuk  or  Savonga-Gambel 
cannot  be  considered  as  reasonable  if  the  major  goal  is  to 
obtain  accurate  estimates  of  the  BST. 

4.  Summary  and  discussion 

In  this  study,  we  comprehensively  analyzed  the  impact 
of  possible  HFR  locations  on  the  accuracy  of  monitoring 
the  MHS  transports  through  the  Bering  Strait  and  their 
further  distribution  in  the  Chukchi  Sea.  The  analysis  is 
based  on  rigorous  estimation  of  the  a  posteriori  error  re¬ 
duction  of  the  transports  through  eight  sections  under  the 
dynamical  constraints  imposed  by  an  OGCM.  This  is 
achieved  using  the  so-called  A-optimality  criterion,  which 
is  based  on  reduction  of  the  trace  of  the  error  covariance. 
We  considered  36  realistic  combinations  of  HFR  pair  de¬ 
ployment  at  existing  Alaskan  settlements. 

Results  of  the  study  indicate  that  a  pair  of  HFRs 
deployed  at  Little  Diomede  Island  and  Cape  Prince  of 
Wales  provides  the  largest  reduction  (fq  =  6.3%  on  aver¬ 
age)  of  the  total  MHS  transport  error.  A  similar  reduction 
(6.1%)  is  obtained  for  HFRs  at  CPW  and  Sinuk. 
Northward-looking  HFR  antennas  in  CPW  and  Little 
Diomede  Island  showed  a  considerably  (10%-20%)  smaller 
error  reduction  in  comparison  to  southward-looking  an¬ 
tennas  at  those  same  sites.  However,  collocation  of  both 
a  northward-  and  southward-looking  antenna  at  DI  or  CPW 
(with  =  5.8%)  may  be  preferable  for  logistic  reasons. 

The  velocity  sensitivity  maps  obtained  in  the  course  of 
the  analysis  indicate  that  placement  of  a  single  mooring 
in  the  eastern  channel  of  the  Bering  Strait  is  best  in 
terms  of  error  reduction  of  the  MHS  transports.  This 
can  be  explained  by  the  fact  that  the  eastern  channel  is 
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much  deeper  than  the  western  one  (Fig.  1)  and  is 
characterized  by  stronger  cross-strait  variations  of  the 
temperature  and  salinity  fields  during  the  ice-free  pe¬ 
riod  considered  here. 

We  also  compared  the  efficiency  of  monitoring  the 
target  transports  by  the  three  moorings  permanently 
maintained  by  NO  A  A  in  the  region  with  the  efficiency 
provided  by  the  optimally  placed  pair  of  HFRs.  Com¬ 
putations  have  shown  a  significant  advantage  of  the 
HFR  monitoring  technique;  the  respective  values  of  fq 
were  6-8  times  larger  for  the  HFRs.  On  the  other  hand, 
moored  observational  platforms  provide  information  in 
the  vertical,  and,  more  importantly,  are  capable  of  op¬ 
eration  under  ice  conditions  when  HFRs  yield  no  mea¬ 
surements.  For  these  reasons,  it  is  highly  desirable  to 
complement  the  existing  moorings  by  HFRs  that  have 
the  potential  to  significantly  increase  the  accuracy  of  TQ 
estimates  during  the  ice-free  season. 

The  technique  used  here  can  also  be  applied  to  opti¬ 
mization  of  multiple  radars  combined  with  permanent 
mooring  deployments.  In  this  case,  the  amount  of  com¬ 
putations  will  increase  significantly,  since  a  mooring  may 
be  placed  at  a  much  larger  number  of  points  than  a  radar 
station.  However,  these  computations  are  quite  feasible 
given  that  they  are  easily  parallelizable  and  inexpensive; 
computation  of  for  a  single  TQ  and  a  given  pair  of 
HFRs  consumed  2.5  CPU  minutes  on  a  laptop. 

In  that  respect,  an  important  issue  is  development  of 
a  general  algorithm  maximizing  gq  with  respect  to  varia¬ 
tions  of  the  observation  operator  D,  whose  structure 
depends  not  only  on  the  location  of  the  observational 
platform  but  on  its  configuration  (such  as  tunable  antenna 
parameters  in  the  case  of  HFRs).  In  the  present  study  we 
immensely  simplified  the  task  by  considering  only  a  lim¬ 
ited  number  of  logistically  feasible  locations  that  allowed 
for  a  solution  to  the  problem  by  brute  force.  In  the  more 
general  formulation,  it  would  be  necessary  to  develop 
effective  algorithms  to  seek  a  maximum  of  fq  with  respect 
to  continuously  varying  parameters. 

Nevertheless,  we  believe  that  the  results  of  this  study 
may  serve  as  guidance  for  future  planning  of  an  HFR 
network  in  the  Bering  Strait  region,  although  permanent 
HFR  positions  also  depend  on  seasonal  ice  cover  and 
other  considerations  not  taken  into  account  in  the  analy¬ 
sis.  It  may  also  be  argued  that  the  results  could  be  model 
dependent.  We  believe,  however,  that  such  dependence 
may  introduce  only  minor  corrections  because  the  July- 
September  reference  state  used  in  our  model  is  robust  and 
statistically  consistent  with  observations. 

In  general,  the  proposed  technique  could  be  applied  to 
the  design  of  observing  systems  in  the  hot  spots  of  the 
World  Ocean  with  typical  scales  (100-200  km)  limited  by 
areal  coverage  of  the  long-range  HFR  systems.  We  believe 


that  increasing  the  efficiency  of  the  observing  systems  in 
these  regions  will  remain  one  of  the  major  scientific  tasks 
in  the  near  future  that  will  promote  further  development 
and  application  of  the  methods  described  in  this  study. 
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